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Abstract 

We review how (dimensionally regulated) scattering amplitudes in M = 4 super- 
Yang-Mills theory provide a useful testing ground for perturbative QCD calcula- 
tions relevant to collider physics, as well as another avenue for investigating the 
AdS/CFT correspondence. We describe the iterative relation for two-loop scatter- 
ing amplitudes in N = 4 super- Yang-Mills theory found in C. Anastasiou et al., 
Phys. Rev. Lett. 91:251602 (2003), and discuss recent progress toward extending it 
to three loops. 



1 Introduction and Collider Physics Motivation 



Maximally supersymmetric (N = 4) Yang-Mills theory (MSYM) is unique 
in many ways. Its properties are uniquely specified by the gauge group, say 
SU(iV c ), and the value of the gauge coupling g. It is conformally invariant for 
any value of g. Although gravity is not present in its usual formulation, MSYM 
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is connected to gravity and string theory through the AdS/CFT correspon- 
dence [1] . Because this correspondence is a weak-strong coupling duality, it is 
difficult to verify quantitatively for general observables. On the other hand, 
such checks are possible and have been remarkably successful for quantities 
protected by supersymmetry such as BPS operators [2], or when an additional 
expansion parameter is available, such as the number of fields in sequences of 
composite, large i?-charge operators [3,4,5,6,7,8]. 

It is interesting to study even more observables in perturbative MSYM, in 
order to see how the simplicity of the strong coupling limit is reflected in the 
structure of the weak coupling expansion. The strong coupling limit should be 
even simpler when the large- N c limit is taken simultaneously, as it corresponds 
to a weakly-coupled supergravity theory in a background with a large radius 
of curvature. There are different ways to study perturbative MSYM. One 
approach is via computation of the anomalous dimensions of composite, gauge 
invariant operators [1,3,4,5,6,7,8]. Another possibility [9], discussed here, is to 
study the scattering amplitudes for (regulated) plane-wave elementary field 
excitations such as gluons and gluinos. 

One of the virtues of the latter approach is that perturbative MSYM scat- 
tering amplitudes share many qualitative properties with QCD amplitudes in 
the regime probed at high-energy colliders. Yet the results and the computa- 
tions (when organized in the right way) are typically significantly simpler. In 
this way, MSYM serves as a testing ground for many aspects of perturbative 
QCD. MSYM loop amplitudes can be considered as components of QCD loop 
amplitudes. Depending on one's point of view, they can be considered either 
"the simplest pieces" (in terms of the rank of the loop momentum tensors in 
the numerator of the amplitude) [10,11], or "the most complicated pieces" in 
terms of the degree of transcendentality (see section 6) of the special functions 
entering the final results [12]. As discussed in section 6, the latter interpreta- 
tion links recent three- loop anomalous dimension results in QCD [13] to those 
in the spin-chain approach to MSYM [5]. 

The most direct experimental probes of short-distance physics are collider 
experiments at the energy frontier. For the next decade, that frontier is at 
hadron colliders — Run II of the Fermilab Tevatron now, followed by startup 
of the CERN Large Hadron Collider in 2007. New physics at colliders always 
contends with Standard Model backgrounds. At hadron colliders, all physics 
processes — signals and backgrounds — are inherently QCD processes. Hence 
it is important to be able to predict them theoretically as precisely as possi- 
ble. The cross section for a "hard," or short-distance-dominated processes, can 
be factorized [14] into a partonic cross section, which can be computed order 
by order in perturbative QCD, convoluted with nonperturbative but measur- 
able parton distribution functions (pdfs). For example, the cross section for 
producing a pair of jets (plus anything else) in a pp collision is given by 
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where s is the squared center-of-mass energy, £ lj2 are the longitudinal (light- 
cone) fractions of the p, p momentum carried by partons a, b, which may be 
quarks, anti-quarks or gluons. The experimental definition of a jet is an in- 
volved one which need not concern us here. The pdf f a (x,fi F ) gives the prob- 
ability for finding parton a with momentum fraction x inside the proton; 
similarly ft, is the probability for finding parton b in the antiproton. The pdfs 
depend logarithmically on the factorization scale fi F , or transverse resolution 
with which the proton is examined. The Mellin moments of f a {x,fiF) are for- 
ward matrix elements of leading-twist operators in the proton, renormalized 
at the scale //p. The quark distribution function q(x,fi), for example, obeys 
/„ dxx j q(x, //) = (p\ [qi + diq] {fl) \p). 



2 Ingredients for a NNLO Calculation 

Many hadron collider measurements can benefit from predictions that are 
accurate to next-to-next-to-leading order (NNLO) in QCD. Three separate 
ingredients enter such an NNLO computation; only the third depends on the 
process: 

(1) The experimental value of the QCD coupling a s {jiRj must be determined 
at one value of the renormalization scale ji R (for example mz), and its 
evolution in /ir computed using the 3-loop /^-function, which has been 
known since 1980 [15]. 

(2) The experimental values for the pdfs / (x, //p) must be determined, ide- 
ally using predictions at the NNLO level, as are available for deep-inelastic 
scattering [16] and more recently Drell-Yan production [17]. The evolu- 
tion of pdfs in /jl f to NNLO accuracy has very recently been completed, 
after a multi-year effort by Moch, Vermaseren and Vogt [13] (previously, 
approximations to the NNLO kernel were available [18]). 

(3) The NNLO terms in the expansion of the partonic cross sections must be 
computed for the hadronic process in question. For example, the parton 
cross sections for jet production has the expansion, 



The quantities A and B have been known for over a decade [19], but C 
has not yet been computed. 



Vab^jjx = a 2 s (A + a s B + a\C + ...). 



(2) 
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Figure 1. LHC Z production [22]. Figure 2. Purely gluonic contributions 

to agg^jjx at NNLO. 



Indeed, the NNLO terms are unknown for all but a handful of collider pro- 
cesses. Computing a wide range of processes at NNLO is the goal of a large 
amount of recent effort in perturbative QCD [20]. As an example of the im- 
proved precision that could result from this program, consider the production 
of a virtual photon, W or Z boson via the Drell-Yan process at the Tevatron 
or LHC. The total cross section for this process was first computed at NNLO 
in 1991 [21]. Last year, the rapidity distribution of the vector boson also be- 
came available at this order [17,22], as shown in fig. 1. The rapidity is defined 
in terms of the energy E and longitudinal momentum p z of the vector boson 
in the center-of-mass frame, Y = A log ( ^Vp* ) ■ ft determines where the vector 
boson decays within the detector, or outside its acceptance. The rapidity is 
sensitive to the x values of the incoming partons. At leading order in QCD, 
X\ = e Y m v /yfs, %2 = e~ Y mv/ y/s, where my is the vector boson mass. 

The LHC will produce roughly 100 million Ws and 10 million Zs per year 
in detectable (leptonic) decay modes. LHC experiments will be able to map 
out the curve in fig. 1 with exquisite precision, and use it to constrain the 
parton distributions — in the same detectors that are being used to search for 
new physics in other channels, often with similar qq initial states. By taking 
ratios of the other processes to the "calibration" processes of single W and Z 
production, many experimental uncertainties, including those associated with 
the initial state parton distributions, drop out. Thus fig. 1 plays a role as a 
"partonic luminosity monitor" [23]. To get the full benefit of the remarkable 
experimental precision, though, the theory uncertainty must approach the 
1% level. As seen from the uncertainty bands in the figure, this precision is 
only achievable at NNLO. The bands are estimated by varying the arbitrary 
renormalization and factorization scales and /ip (set to a common value 
/i) from my i "2 to 2my. A computation to all orders in a s would have no 
dependence on /i. Hence the //-dependence of a fixed order computation is 
related to the size of the missing higher-order terms in the series. Although 
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sub-1% uncertainties may be special to W and Z production at the LHC, 
similar qualitative improvements in precision will be achieved for many other 
processes, such as di-jet production, as the NNLO terms are completed. 



Even within the NNLO terms in the partonic cross section, there are several 
types of ingredients. This feature is illustrated in fig. 2 for the purely gluonic 
contributions to di-jet production, a gg ^jjx- In the figure, individual Feynman 
graphs stand for full amplitudes interfered ( x ) with other amplitudes, in order 
to produce contributions to a cross section. There may be 2, 3, or 4 partons in 
the final state. Just as in QED it is impossible to define an outgoing electron 
with no accompanying cloud of soft photons, also in QCD sensible observables 
require sums over final states with different numbers of partons. Jets, for 
example, are defined by a certain amount of energy into a certain conical 
region. At leading order, that energy typically comes from a single parton, 
but at NLO there may be two partons, and at NNLO three partons, within 
the jet cone. 



Each line in fig. 2 results in a cross-section contribution containing severe 
infrared divergences, which are traditionally regulated by dimensional regula- 
tion with D = 4 — 2e. Note that this regulation breaks the classical conformal 
invariance of QCD, and the classical and quantum conformal invariance of 
M = 4 super- Yang-Mills theory. Each contribution contains poles in e ranging 
from 1/e 4 to 1/e. The poles in the real contributions come from regions of 
phase-space where the emitted gluons are soft and/or collinear. The poles in 
the virtual contributions come from similar regions of virtual loop integra- 
tion. The virtual x real contribution obviously has a mixture of the two. The 
Kinoshita-Lee-Nauenberg theorem [24] guarantees that the poles all cancel 
in the sum, for properly-defined, short- distance observables, after renormal- 
izing the coupling constant and removing initial-state collinear singularities 
associated with renormalization of the pdfs. 



A critical ingredient in any NNLO prediction is the set of two-loop ampli- 
tudes, which enter the doubly- virtual x real interference in fig. 2. Such ampli- 
tudes require dimensionally-regulated all-massless two-loop integrals depend- 
ing on at least one dimensionless ratio, which were only computed beginning in 
1999 [25,26,27]. They also receive contributions from many Feynman diagrams, 
with lots of gauge-dependent cancellations between them. It is of interest to 
develop more efficient, manifestly gauge-invariant methods for combining di- 
agrams, such as the unitarity or cut-based method successfully applied at one 
loop [10] and in the initial two- loop computations [28]. 
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Figure 3. Illustration of soft-collinear (left) and pure-collinear (right) one-loop di- 
vergences. 

3 Af = 4 Super- Yang-Mills Theory as a Testing Ground for QCD 



M = 4 super- Yang-Mills theory serves an excellent testing ground for pertur- 
bative QCD methods. For n-gluon scattering at tree level, the two theories 
in fact give identical predictions. (The extra fermions and scalars of MSYM 
can only be produced in pairs; hence they only appear in an n-gluon ampli- 
tude at loop level.) Therefore any consequence of M = 4 super symmetry, such 
as Ward identities among scattering amplitudes [29], automatically applies 
to tree-level gluonic scattering in QCD [30]. Similarly, at tree level Witten's 
topological string [31] produces MSYM, but implies twistor-space localization 
properties for QCD tree amplitudes. (Amplitudes with quarks can be related 
to supersymmetric amplitudes with gluinos using simple color manipulations.) 



3.1 Pole Structure at One and Two Loops 



At the loop- level, MSYM becomes progressively more removed from QCD. 
However, it can still illuminate general properties of scattering amplitudes, in 
a calculationally simpler arena. Consider the infrared singularities of one-loop 
massless gauge theory amplitudes. In dimensional regularization, the leading 
singularity is 1/e 2 , arising from virtual gluons which are both soft and collinear 
with respect to a second gluon or another massless particle. It can be char- 
acterized by attaching a gluon to any pair of external legs of the tree-level 
amplitude, as in the left graph in fig. 3. Up to color factors, this leading diver- 
gence is the same for MSYM and QCD. There are also purely collinear terms 
associated with individual external lines, as shown in the right graph in fig. 3. 
The pure-collinear terms have a simpler form than the soft terms, because 
there is less tangling of color indices, but they do differ from theory to theory. 

The full result for one-loop divergences can be expressed as an operator 
which acts on the color indices of the tree amplitude [32] . Treating the L-loop 
amplitude as a vector in color space, \A^), the one-loop result is 

\A^) = I^(e)\A^} + \A^}, (3) 
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where |^4^ ,fin ) is finite as e — > 0, and 
I«(e) 



on 



1 e c 

2fO^)gg Ti - Tj 



e 2 + T? 6 



j4_ 



(4) 



where 7 is Euler's constant and Sjj = (fcj + fcj) 2 is a Mandelstam invariant. 
The color operator Tj • T 3 - = T^T? and factor of (/^/j/ ( — s ij)) e arise from soft 
gluons exchanged between legs i and j, as in the left graph in fig. 3. The pure 
1/e poles terms proportional to 7, have been written in a symmetric fashion, 
which slightly obscures the fact that the color structure is actually simpler. 
We can use the equation which represents color conservation in the color-space 
notation, Y^=i Tj = 0, to simplify the result. At order 1/e we may neglect the 
— s ij)Y factor in the 7, terms, and we have Z^Tj • Tj %/T 2 = —7$. 
So the color structure of the pure 1/e term is actually trivial. For an n-gluon 
amplitude, the factor 7$ is set equal to its value for gluons, which turns out to 
be 7 9 = 60, the one-loop coefficient in the /3-function. Hence the pure-collinear 
contribution vanishes for MSYM, but not for QCD. 

The divergences of two-loop amplitudes can be described in the same for- 
malism [32]. The relation to soft-collinear factorization has been made more 
transparent by Sterman and Tejeda-Yeomans, who also predicted the three- 
loop behavior [33]. Decompose the two-loop amplitude as 

l-4 2) > = I (2) (c)K 0) > +I {1 \e)\AU) + \A^) , (5) 
where \A^' &n ) is finite as e — > and 

I<%) = (l« M + ^) + e "';[ 1 (1 _- 2t) + /<) *•><*) 

„e7 r n n fj( 2 ) , ,2 n i 

+ * . - V E T, • T.-^V f — ) + H (2) • (6) 
4er(l-e)[ ^ J T 2 U^J J ^ ' 

(2) (2) 

Here -K" and ifj depend on the theory, and H\ , like 7,, also depends on 
the external leg i For QCD, if has long been known from soft-gluon resum- 
mation [34], while were found by explicit computation of four-parton 
two-loop scattering amplitudes [35,36,37]. For MSYM, the quantities are nat- 
urally simpler, 



K 



-C2 c A , 



H (2),Af=4 = (3^,2 



(7) 
(8) 
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where Ca = N c is the adjoint Casimir value. The quantity IT 2 ) has non- 
trivial, but purely subleading-in-iV c , color structure. It is associated with soft, 
rather than collinear, momenta [37,33], so it is theory-independent, up to color 
factors. An ansatz for it for general n has been presented recently [38]. 



3.2 Recycling Cuts in MSYM 



An efficient way to compute loop amplitudes, particularly in theories with 
a great deal of supersymmetry, is to use unitarity and reconstruct the am- 
plitude from its cuts [10,38]. For the four-gluon amplitude in MSYM, the 
two-loop structure, and much of the higher-loop structure, follows from a sim- 
ple property of the one-loop two-particle cut in this theory. For simplicity, we 
strip the color indices off of the four-point amplitude Af\ by decomposing 
it into color-ordered amplitudes Af\ whose coefficients are traces of SU(iV c ) 
generator matrices (Chan-Paton factors), 



Af\k 1 ,a 1 ;k 2 ,a 2 ;k 3 ,a 3 ;h,a 4 )=g 2 ]T Ti (T a ^T a ^T a ^T a ^) 

X^4 0) Op(l) ? fcp(2), fcp(3), fcp(4)) • (9) 

The two-particle cut can be written as a product of two four-point color- 
ordered amplitudes, summed over the pair of intermediate Af = 4 states S, S' 
crossing the cut, which evaluates to 



]T Af\k x , k 2 ,£ s , —t's') x Af\£> s „-l s ,k 3 ,k 4 ) 

S,S'eAT=4 

= *s 12 s 23 Af ) {k 1 , k 2 , k 3 , h) x - , (10) 

where £' = £ — k\ — k 2 . This equation is also shown in fig. 4. The scalar 
propagator factors in eq. (10) are depicted as solid vertical lines in the figure. 
The dashed line indicates the cut. Thus the cut reduces to the cut of a scalar 
box integral, defined by 

r d 4 ~ 2t f 1 
4 J {27iy-^£ 2 {£-k l ) 2 {£-k 1 -k 2 ) 2 {£ + k A y y } 



One of the virtues of eq. (10) is that it is valid for arbitrary external states 
in the M = 4 multiplet, although only external gluons are shown in fig. 4. 
Therefore it can be re-used at higher loop order, for example by attaching yet 
another tree to the left. 
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Figure 4. The one-loop two-particle cuts for the four-point amplitude in MSYM 
reduce to the tree amplitude multiplied by a cut scalar box integral (for any set of 
four external states). 



N=4 



l S 12 S 2 3 




+ perms 



Figure 5. The two-loop gg — > gg amplitude in MSYM [11,39]. The blob on the 
right represents the color-ordered tree amplitude A± . (The quantity si2S23^4 
transforms symmetrically under gluon interchange.) In the the brackets, black lines 
are kinematic 1/p 2 propagators, with scalar (4> 3 ) vertices. Green lines are color 8 ab 
propagators, with structure constant (f abc ) vertices. The permutation sum is over 
the three cyclic permutations of legs 2,3,4, and makes the amplitude Bose symmetric. 



At two loops, the simplicity of eq. (10) made it possible to compute the 
two-loop gg — > gg scattering amplitude in that theory (in terms of specific 
loop integrals) in 1997 [11], four years before the analogous computations in 
QCD [36,37]. All of the loop momenta in the numerators of the Feynman di- 
agrams can be factored out, and only two independent loop integrals appear, 
the planar and nonplanar scalar double box integrals. The result can be writ- 
ten in an appealing diagrammatic form, fig. 5, where the color algebra has the 
same form as the kinematics of the loop integrals [39]. 

At higher loops, eq. (10) leads to a "rung rule" [11] for generating a class 
of (L + l)-loop contributions from L-loop contributions. The rule states that 
one can insert into a L-loop contribution a rung, i.e. a scalar propagator, 
transverse to two parallel lines carrying momentum t\ + £2, along with a factor 
of i(£i + £2) 2 m the numerator, as shown in fig. 6. Using this rule, one can 
construct recursively the external and loop-momentum-containing numerators 
factors associated with every </> 3 -type diagram that can be reduced to trees by a 
sequence of two-particle cuts, such as the diagram in fig. 7 a. Such diagrams can 
be termed "iterated 2-particle cut-constructible," although a more compact 
notation might be 'Mondrian' diagrams, given their resemblance to Mondrian's 
paintings. Not all diagrams can be computed in this way. The diagram in fig. 7b 
is not in the 'Mondrian' class, so it cannot be determined from two-particle 
cuts. Instead, evaluation of the three-particle cuts shows that it appears with 
a non-vanishing coefficient in the subleading-color contributions to the three- 
loop MSYM amplitude. 
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Figure 6. The rung rule for MSYM. Figure 7. (a) Example of a 'Mondrian' 

diagram which can be determined re- 
cursively from the rung rule, (b) The 
first non-vanishing, non-Mondrian dia- 
grams appear at three loops in nonplanar, 
subleading-color contributions. 

4 Iterative Relation in M = 4 Super- Yang-Mills Theory 



Although the two-loop gg — > gg amplitude in MSYM was expressed in terms 
of scalar integrals in 1997 [11], and the integrals themselves were computed 
as a Laurent expansion about D = 4 in 1999 [25,26], the expansion of the 
M = 4 amplitude was not inspected until last fall [9], considerably after similar 
investigations for QCD and Af = 1 super- Yang-Mills theory [36,37]. It was 
found to have a quite interesting "iterative" relation, when expressed in terms 
of the one-loop amplitude and its square. 

At leading color, the L-loop gg — ► gg amplitude has the same single-trace color 
decomposition as the tree amplitude, eq. (9). Let AQ be the ratio of this 
leading-color, color-ordered amplitude to the corresponding tree amplitude, 
omitting also several conventional factors, 



(L),Af=4 planar 



2e~^g 2 N c 



(4tt) 



L 



v# x M[ L) . (12) 



Then the iterative relation (see also fig. 8) is 

Mf(e) = ^(M 4 (1) (e)) 2 + /(e) Mf(2e) - V 2 ) 2 + 0(e) , (13) 



where /(e) = (^(1 - e) - ^(l))/e = + C 3 e + Qe 2 + ■■•)• 

The analogous relation for gluon-gluon scattering at two loops in QCD takes 
a similar form at the level of the pole terms in e, thanks to the general re- 
sult (5). But the finite remainder —\{C,2) 2 is replaced by approximately six 
pages of formulas (!), including a plethora of polylogarithms, logarithms and 
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N=4 planar 





KC2) 2 + O(e) 




(C2 + <3 + e 2 C4 + ...) 



Figure 8. Schematic depiction of the iterative relation (13) between two-loop and 
one-loop MSYM amplitudes. 

polynomials in ratios of invariants s/t, s/u and t/u [37]. The polylogarithm is 
defined by 



It appears with degree m up to 4 at the finite, order e°, level; and up to 
degree 4 — i in the 0(e~ l ) terms. In the case of MSYM, identities relating 
these polylogarithms are needed to establish eq. (13). 

Although the O(e ) term in eq. (13) is miraculously simple, as noted above 
the behavior of the pole terms is not a miracle. It is dictated in general terms 
by the cancellation of infrared divergences between virtual corrections and 
real emission [24]. Roughly speaking, for this cancellation to take place, the 
virtual terms must resemble lower-loop amplitudes, and the real terms must 
resemble lower-point amplitudes, in the soft and collinear regions of loop or 
phase-space integration. 

At the level of the finite terms, the iterative relation (13) can be understood 
in the Regge/BFKL limit where s ^> t, because it then corresponds to expo- 
nentiation of large logarithms of s/t [40]. For general values of s/t, however, 
there is no such argument. 

The relation is special to D = 4, where the theory is conformally invariant. 
That is, the 0(e x ) remainder terms cannot be simplified significantly. For ex- 
ample, the two-loop amplitude M4 (e) contains at 0(e 1 ) all three independent 
Li 5 functions, Li 5 (— s/u), Li 5 (— t/u) and Li 5 (— s/t), yet [M4 (e)] 2 has only the 
first two of these [9]. 

The relation is also special to the planar, leading-color limit. The subleading 
color-components of the finite remainder l^ 2 - 1 ' 611 ) defined by eq. (5) show no 
significant simplification at all. 




X 



ln(l — x). 



(14) 
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For planar amplitudes in the D — ■> 4 limit, however, there is evidence that 
an identical relation also holds for an arbitrary number n of external legs, 
at least for certain "maximally helicity-violating" (MHV) helicity amplitudes. 
This evidence comes from studying the limits of two-loop amplitudes as two of 
the n gluon momenta become collinear [9,38,41]. (Indeed, it was by analyzing 
these limits that the relation for n = 4 was first uncovered.) The collinear 
limits turn out to be consistent with the same eq. (13) with M4 replaced by 
M n everywhere [9], i.e. 

Mi 2 \e) = \{M^{e)) 2 + /(e) M«(2e) - \{C 2 f + 0(e) . (15) 

The collinear consistency does not constitute a proof of eq. (15), but in light 
of the remarkable properties of MSYM, it would be surprising if it were not 
true in the MHV case. Because the direct computation of two-loop amplitudes 
for n > 4 seems rather difficult, it would be quite interesting to try to examine 
the twistor-space properties of eq. (15), along the lines of refs. [31,42]. (The 
right-hand-side of eq. (15) is not completely specified at order 1/e and e° for 
n > 4. The reason is that the order e and e 2 terms in M^e), which contribute 
to the first term in eq. (15) at order 1/e and e°, contain the D = 6 — 2e 
pentagon integral [43], which is not known in closed form. On the other hand, 
the differential equations this integral satisfies may suffice to test the twistor- 
space behavior. Or one may examine just the finite remainder M^ L )' fin defined 
via eq. (5).) 

It may soon be possible to test whether an iterative relation for planar MSYM 
amplitudes extends to three loops. An ansatz for the three-loop planar gg — > gg 
amplitude, shown in fig. 9, was provided at the same time as the two-loop re- 
sult, in 1997 [11]. The ansatz is based on the "rung-rule" evaluation of the 
iterated 2-particle cuts, plus the 3-particle cuts with intermediate states in 
D = 4; the 4-particle cuts have not yet been verified. Two integrals, each be- 
ginning at (9(e~ 6 ), are required to evaluate the ansatz in a Laurent expansion 
about D = 4. (The other two integrals are related by s <-> t.) The triple ladder 
integral on the top line of fig. 9 was evaluated last year by Smirnov, all the 
way through O(e ) [44]. Evaluation of the remaining integral, which contains 
a factor of (£ + k^) 2 in the numerator, is in progress [45]; all the terms through 
0(e~ 2 ) agree with predictions [33], up to a couple of minor corrections. 



5 Significance of Iterative Behavior? 

It is not yet entirely clear why the two-loop four-point amplitude, and prob- 
ably also the n-point amplitudes, have the iterative structure (15). However, 
one can speculate that it is from the need for the perturbative series to 



12 




Figure 9. Graphical representation of the three-loop amplitude for MSYM in the 
planar limit. 

be summable into something which becomes "simple" in the planar strong- 
coupling limit, since that corresponds, via AdS/CFT, to a weakly-coupled 
supergravity theory. The fact that the relation is special to the conformal 
limit D — > 4, and to the planar limit, backs up this speculation. Obviously 
it would be nice to have some more information at three loops. There have 
been other hints of an iterative structure in the four-point correlation func- 
tions of chiral primary (BPS) composite operators [46], but here also the exact 
structure is not yet clear. Integrability has played a key role in recent higher- 
loop computations of non-BPS spin-chain anomalous dimensions [4,5,6,8]. By 
imposing regularity of the BMN 'continuum' limit [3], a piece of the anoma- 
lous dimension matrix has even been summed to all orders in g 2 N c in terms 
of hypergeometric functions [7] . The quantities we considered here — gauge- 
invariant, but dimensionally regularized, scattering amplitudes of color non- 
singlet states — are quite different from the composite color-singlet operators 
usually treated. Yet there should be some underlying connection between the 
different perturbative series. 



6 Aside: Anomalous Dimensions in QCD and MSYM 



As mentioned previously, the set of anomalous dimensions for leading-twist 
operators was recently computed at NNLO in QCD, as the culmination of a 
multi-year effort [13] which is central to performing precise computations of 
hadron collider cross sections. Shortly after the Moch, Vermaseren and Vogt 
computation, the anomalous dimensions in MSYM were extracted from this 
result by Kotikov, Lipatov, Onishchenko and Velizhanin [12]. (The MSYM 
anomalous dimensions are universal; supersymmetry implies that there is only 
one independent one for each Mellin moment j.) This extraction was non- 
trivial, because MSYM contains scalars, interacting through both gauge and 
Yukawa interactions, whereas QCD does not. However, Kotikov et al. noticed, 
from comparing NLO computations in both leading-twist anomalous dimen- 
sions and BFKL evolution, that the "most complicated terms" in the QCD 
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computation always coincide with the MSYM result, once the gauge group 
representation of the fermions is shifted from the fundamental to the adjoint 
representation. One can define the "most complicated terms" in the x-space 
representation of the anomalous dimensions — i.e. the splitting kernels - 
as follows: Assign a logarithm or factor of n a transcendentality of 1, and a 
polylogarithm Li m or factor of ( m = Li m (l) a transcendentality of m. Then 
the most complicated terms are those with leading transcendentality. For the 
NNLO anomalous dimensions, this turns out to be transcendentality 4. (This 
rule for extracting the MSYM terms from QCD has also been found to hold 
directly at NNLO, for the doubly-virtual contributions [38].) Strikingly, the 
NNLO MSYM anomalous dimension obtained for j = 4 by this procedure 
agrees with a previous result derived by assuming an integrable structure for 
the planar three-loop contribution to the dilatation operator [5]. 



7 Conclusions and Outlook 



M = 4 super- Yang-Mills theory is an excellent testing ground for techniques 
for computing, and understanding the structure of, QCD scattering amplitudes 
which are needed for precise theoretical predictions at high-energy colliders. 
One can even learn something about the structure of M = 4 super- Yang-Mills 
theory in the process, although clearly there is much more to be understood. 
Some open questions include: Is there any AdS/CFT "dictionary" for color 
non-singlet states, like plane- wave gluons? Can one recover composite operator 
correlation functions from any limits of multi-point scattering amplitudes? Is 
there a better way to infrared regulate M = 4 supersymmetric scattering 
amplitudes, that might be more convenient for approaching the AdS/CFT 
correspondence, such as compactification on a three-sphere, use of twistor- 
space, or use of coherent external states? Further investigations of this arena 
will surely be fruitful. 
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